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Abstract 



The possibilities to determine the flavor structure of the polarized sea (antiquark) 
distributions of the nucleon via vector boson (7*, W , Z°) production at high 
energy polarized hadron-hadron colliders, such as the Relativistic Heavy-Ion Col- 
lider (RHIC), are studied in detail. In particular the perturbative stability of the 
expected asymmetries in two representative models for the (un)broken flavor struc- 
ture are investigated by confronting perturbative QCD leading order predictions of 
the expected asymmetries with their next-to-leading order counterparts. 



1 Introduction 



The determination of the polarized parton content of the nucleon via measurements of 
the inclusive structure functions gi' n {x, Q 2 ) does not provide detailed information con- 
cerning the flavor structure of these distributions, in complete analogy to unpolarized 
deep inelastic structure functions. In particular the flavor structure of the antiquark (sea) 
distributions is not fixed and one needs to resort to semi-inclusive deep inelastic hadron 
production for this purpose |], [J. The resulting antiquark distributions Aq are, however, 
not reliably determined by this method @, |J for the time being due to their dependence 
on the rather poorly known quark fragmentation functions at low scales. 

A more reliable determination is provided via inclusive vector boson (7*, W ± , Z°) pro- 
duction in polarized hadron-hadron collisions as envisaged at RHIC (BNL) or at HERA-N 
(DESY) whose potential for the determination of Am and Ad will be studied in this paper 
in leading (LO) as well as next-to-leading (NLO) order of QCD. 

Recently the LO Drell-Yan (7*) dilepton production process for future polarized pp 
and pd collision experiments has been suggested and studied |5], || for delineating the 
flavor-asymmetry of the polarized light sea distributions Au(x, Q 2 ) 7^ Ad(x,Q 2 ). In 
particular, the relativistic field theoretic chiral quark-soliton model JF], || |J as well as 
a more recent analysis based on the statistical parton model [[l(J predict \Ad(x, Q 2 )\ > 
Au(x,Q 2 ). Similar expectations [[TTJ, |7|, [T2L (TO, |TB|, partly from first principles, hold also 



for the unpolarized light sea distributions, d > u, which have been already confirmed, 
as is well known, by Drell-Yan production |14| pp and pd experiments | 15[j . In 



general, the flavor-asymmetry of the light sea distributions can be understood in terms of 
flavor mass asymmetries and 'Pauli-blocking' effects being related to the Pauli exclusion 



principle |16|, O, [13 



In order to analyze the determination of (flavor-asymmetric) light antiquark (sea) 
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distributions one needs some alternative models for Aq, introduced in Section 2, and 
the possibility of their experimental distinction will be investigated in Section 3. Our 
conclusions are drawn in Sect. 4 and the Appendix contains all expressions of LO and 
NLO cross sections relevant for our calculations. 



2 Models for the Polarized Antiquark Distributions 
of the Nucleon 



(i) Standard unbroken sea scenario: Here one assumes, as in most analyses of polarization 
data performed thus far, a flavor symmetric sea, i.e. 

Au(x, Q 2 ) = Ad(x, Q 2 ) = As(x, Q 2 ) = Aq(x, Q 2 ) (2.1) 

where as usual Au sea = Au, Ad sca = Ad and As = As. The adopted LO and NLO 



in 



distributions will be taken from the recent analysis of the AA Collaboration JT8 
particular their LO and NLO-1 ones. 

(ii) Broken sea scenario: Here one assumes Au(x,Q 2 ) ^ Ad(x,Q 2 ) ^ As(x,Q 2 ), i.e. a 
broken flavor symmetry as motivated by the situation in the corresponding unpolarized 



sector [15, O]. As mentioned in the Introduction, present polarization data do not pro- 



vide detailed and reliable information concerning flavor symmetry breaking and we shall 
therefore utilize antiquark distributions extracted via the phenomenological ansatz of [pjj 
which is confirmed in the unpolarized sector and which moreover agrees well with the 
predictions obtained within the framework of the chiral quark-soliton model []7|, £| 0] 
and with a recent analysis of semi-inclusive deep inelastic data 0]. 

These flavor-asymmetric distributions, henceforth denoted by Af, are related to the 



(flavor-symmetric) AAC |18| distributions Af by the relations 
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(2.2) 

with Aq v = Aq — Aq and Aq given by (2.1) at the input scale Ql = 1 GeV 2 , and where 

Am - Ad 



-Aq 



Au + Ad - 2Aq 
as follows from the Pauli-blocking relation [|T7|] 

Ad'(x, Ql)/Au'(x, Ql) = Au'(x, Q 2 )/Ad'(x, Ql) 

combined with the constraints (q = u, d) 

Aq\x,Ql)+Aq!{x,Ql) = Aq(x } Ql) + Aq(x } Ql) 



(2.3) 



(2.4) 



^<i>,Ql) = Yl A( iv( x 'Q 

q=u,d q=u,d 



(2.5) 



needed [19 to preserve the quality of the fit |18|] to g^' n (x,Q 2 ) within the standard un- 
broken sea scenario. Note that Au' — Ad' = 20 and that the 'breaking' function <p(x) in 
(2.3) can be simply parametrized in LO and NLO as well, 



x 



0.536 



(0.27 - O.Ux + 13.57x 2 - 10.17x 3 )(l - x) 



10.5 



X(p{X) LO 

x0(x) NLO = a; a32 (0.12 + 0.90x - 6.68a; 2 + 20.06x 3 - 20.36a; 4 )(l - x) 7 - 7 , (2.6) 
which allows to form Mellin n-moments of needed for performing the Q 2 -evolutions in 



n-moment space as described, for example, in ||20|| . The resulting distributions are shown 
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in fig. 1 at the input scale Ql as well as at Q 2 = 25 GeV 2 and M^, and for comparison 



we also show the flavor-symmetric sea density Aq in (2.1) of AAC |18 |. It should be 
emphasized that our broken distributions should not be used well below x ~ 10 -2 since 
our breaking-ansatz |T7|, [19]] Ad'/Au' = (Aq - <j))/(Aq + 4>) = (Au - 4>)/(Ad + 0) which 



results in (2.3) gives rise to artificial oscillations ]nj below x = 0.01 due to the appearance 
of differences of parton distributions in (2.3). All our subsequent analyses employ these 
broken distributions well above x = 10~ 2 . (For obvious reasons we skip from now on the 
'prime' notation in the above equations for the flavor-asymmetric distributions and simply 
refer to 'AAC in connection with the original AAC-densities with their flavor-symmetric 
sea Aq.) 

Furthermore, the unpolarized cross sections needed for the asymmetry calculations of 
Sect. 3 will be evaluated with the LO and NLO unpolarized parton distributions GRV98 



M utilized by the AAC p . 



3 Cross Section Asymmetries for Inclusive Hadronic 
Vector Boson Production 



As stated in the Introduction, the inclusive hadronic vector boson production yields reli- 
able information on the antiquark (sea) content of the nucleon. We shall consider double 
(single) asymmetries as obtained with doubly (singly) polarized hadron beams. The rel- 
evant expressions for the corresponding cross sections in LO and NLO are collected, for 
convenience, in the Appendix. 

We shall successively study 7* and W ± , Z° production in Sections 3.1 and 3.2, respec- 
tively, starting with estimates of expected production rates at RHIC (BNL) and turning 
then to a study of the optimal ways to gain insight about the flavor structure of the 
antiquark distributions in the nucleon. 
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3.1 7* Production 



The expected production rates at RHIC are determined via the following expected energies 
and integrated luminosities |2f: y/S = 50 - 500 GeV, Cpp(y/S) ~ (y/S/500 GeV) 



800 pb~ l . The polarization rate for the proton beams is expected to be P^ ~ 0.7 and will be 
obviously lower for polarized neutrons, i.e. deuteron beams for example. In our estimates 
for the expected statistical errors we shall use P = 0.7 everywhere which represents an 
idealized underestimate of these rather crucial errors for the case of d(n) beams. 

Asymmetries involving only polarized protons, although easiest to be measured ex- 
perimentally, are unlikely to provide unique signatures for a flavor-broken polarized light 

sea: 

lo 4 [Au(l)Au(2) + Au(l)Au(2)] + Ad(l)Ad(2) + Ad(l)Ad(2) 
4 [u(l)u(2) + u(l)u(2)] + d(l)J(2) + d(l)d(2) 

where (A) a denotes the relevant differential cross sections summarized in the Appendix 
and where we make use of the simple abbreviated '1,2' notation introduced in the Ap- 
pendix as well. The size and even the sign of this asymmetry depends already strongly 



on the choice of the particular set of unbroken polarized parton densities [23|. For our 
sets of (un)broken distributions in Sect. 2 we obtain, for example, at = 50 GeV and 
for an invariant dilepton mass M = 5 GeV in pp — > r y*(M)X — > fi + fi~~X an asymmetry 
Alp ~ +2% for the unbroken (AAC) densities and Al^ ~ —(4 to 6)% for the flavor-broken 
scenario for all values of xf- 



As in the case of flavor-broken unpolarized parton distributions [14, 15| , additional 



polarized pn{pd) reactions are required for delineating the flavor structure of the polarized 
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light (anti) quark sea Au and Ad: 



AaZ^ ± AaZ~ 

4±n - -^r— ^ (3-2) 
lo [4Au(l) ± Ad{l)\ [Au{2) ± Ad(2)] + [4Au(l) ± Ad{l)\ [Au{2) ± Ad(2)] 



Xf— >1 



[4u(l) + d(l)] [u{2) + d(2)] + [4m(1) + d(l)} [u(2) + d{2)} 
[4Au(l) ± Ad{l)] [Au{2) ± Ad{2)] 



[Au(l) + d(l)\ [i2(2) + d(2)] 



*W e= i(l±A s Z;/A^) (3.3) 

lo 1 [4A«(1) ± Ad(l)] [Aw(2) ± Ad(2)] + [4Aw(l) ± Ad(l)] [Au(2) ± Ad(2)] 
2 4[Am(1)Au(2) + Au(1)Am(2)] + Ad(l)Ad(2) + Ad(l)Ad(2) 



XF ^i 1 [4Au(l) ± Ad(l)] [Au(2) ± Ad(2)] _ 1 
* 2 4Aw(l)Aw(2) + Ad(l)Ad(2) ~ 2 



M2) 
Au(2) 



where for xp — > 1, i.e. xi — > 1 and £2 — > 0, the small Ag(l) terms can be neglected 
and in the very last 'crude' (LO) approximation of R p ± n the Ad(l) terms are neglected 
with respect to 4Au(l). This latter ratio of polarized cross sections in (3.3) has been 
suggested and studied in LO originally in as a very sensitive observable for testing the 
flavor-asymmetry of the polarized light sea. Indeed, at small x, R p ± n is sensitive directly 
to the ratio Ad/ Au whereas AZ ±ii in (3.2) is proportional to the isoscalar and isovector 
combinations Am ± Ad. 

First we present in fig. 2 the unpolarized cross section dajl/dM dxp for \fS = 50 and 
100 GeV. It is obvious that only for \xp\ £ 0.5 and not too large dilepton masses M useful 
production rates can be obtained. From 

1 



4,2 = 2 [\/ X F + AM2 I S ± XF ) ( 3 - 4 ) 

one infers that the small-X2 region, relevant for our study of sea (antiquark) distributions, 
implies small r = M 2 /S values. Here, furthermore, the NLO contributions are genuine 
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0(a s ) corrections and small. Since the relevant asymmetries decrease with yS we have 
chosen v^S* = 50 GeV and M = 5 - 8 GeV to be the appropriate range for our study of 
asymmetries related to Drell-Yan dilepton pairs notwithstanding the fact, shown in fig. 2, 
that the unpolarized cross sections do increase with \/S. 

As expected, AZ_ H and in particular R pi i = R p+n are the best indicators for the flavor 
structure of the antiquark (sea) distibutions: Fig. 3 presents our LO and NLO results with 
the latter ones being, furthermore, remarkably stable with respect to sizeable variations 
of the factorization scale fiF- It should be noted that for all scenarios of polarized parton 
distributions with a flavor-symmetric light sea, Au = Ad, we have R p d — > 1 and A^_^ — > 
as x p — ^ 1 a limit which is already reached for xp^ 0.2 for most sets of polarized parton 
densities as illustrated for AAC in fig. 3. The statistical errors shown in fig. 3 are obtained 



from (see, e.g., |22], |2J) 



A * p± " K ± ^ TPS? T^b' (3 ' 6) 

assuming Pi = P 2 = P = 0.7 for the beam polarizations, £ = 80 pb^ 1 and bin widths 
Axf = 0.1, AM = 1 GeV for calculating bin-integrated cross sections. 



3.2 W ± and Z° Production 



The production of these vector bosons via (un)polarized pp(n) collisions, pp{n) — > 
W ± X — * /i ± n X and pp{n) — > Z°X — > fi + fi~X, affords of course higher cm. en- 
ergies which are, however, available at RHIC (v^Smax = 500 GeV) |22| . The expected 
cross sections, presented in fig. 4, are comparable to those for 7* production at = 50 
GeV in fig. 2. These cross sections are partly more sensitive to the flavor structure of the 
light sea, although at far larger scales fip ~ Myy± z°, an d may discern its polarized flavor 
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structure not only in doubly polarized collisions but also in the single polarization mode 
[25|-||27|| where only one beam is polarized. Possible benefits of this latter mode are the 



expected lower statistical errors 

AA 



1 



(3.7) 



i.e., a factor P{~ , as compared to (PiP 2 )~ 1 for the doubly polarized mode in (3.5), as 
well as possibly higher luminosities. As seen in fig. 4, the relevant rapidity range is \y\£> 1 
which covers the interesting range of Bjorken-x for the sea distributions, 0.06 £ x <, 0.4 . 



The relevant single and double helicity asymmetries are 



A 



VP 



aW- 
A PP 



A 



w+ 

VP 



= Aa^ /a w 

lo -A«(1)J(2) + Ad(l)«(2) 
w(l)d(2) + J(l)u(2) 

= A(T W ~ I (T W ~ 

pp I pp 

lo -Ad(l)u(2) + A«(l)d(2) 
d(l)u(2) +u(l)d(2) 

= Aa^- + /^ + 
— pp I pp 

Au(l)Ad(2) + Ad(l)Au(2) 



-A«(l)/«(1), y>+± 

AJ(1)M1), y<-§ 



-Ad(l)/d(l), y>+| 
A«(l)/«(1), y<-§ 



LO 



u(l)d(2) + d(l)u(2) 



Au(l) Ad(2) 
' u(l) J(2) 



.4 



w- 



/>/< 



LO 



Ad(l)A«(2) + A«(l)AtZ(2) j/>+± Ad(l) Au(2) 

* d(l) u(2) 



d(l)ii(2) +u(l)d(2) 



(3- 



(3.9) 



(3.10) 



(3.11) 



with obvious generalizations to AWj and A^ , and where 



±v 



(3.12) 
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In addition we also study the following ratios of singly and doubly polarized cross sections 



Aaf p + lo -A«(l)d(2)+Ad(l)w(2) y<-| Ad(l) u{2) 



Act™- -Ad(l)u(2) + Au(l)d(2) A«(l) d(2) 



(3.13) 



Lg -A^(l)[^(2) +J(2)] + AJ(1)K2) +d(2)] y<-j Agl) 
A(7^" -Ad(l)[w(2) + d(2)] + A«(l)K2) + d(2)j * Au(l) 1 ' J 



Aa^- + lo Au(l)Ad(2) + Ad(l)A«(2) y>+§ Au(l) Ad(2) 

ActKT ~ Ad(l)Au(2) + Au(l)Ad(2) ' Ad(l) Au(2) 



(3.15) 



Aa ^- + Lg M1MM2) + Ad(2)] + Ad(l)[Att(2) + rf(2)] y <-^ ; Agl) 

Ao-^-T ~ Ad(l)[Au(2) + Ad(2)] + Aw(l)[Au(2) + Ad(2)] * Au(l) ' [ ' ' 

The singly polarized asymmetries in (3.8) and (3.9) are, for negative values of y, dominated 
by the polarized antiquark distributions. (Note that y — — 1 corresponds to — 0.06.) 
The limiting values Aq/q and Aq/ q of at the scale n 2 F = are shown in fig. 5 where 
they are compared with the flavor-unbroken (AAC) scenario as well. It is conceivable that 
such differences can be delineated by future RHIC experiments taken into account their 
expected statistical accuracy p2| . 



In fig. 6 we present the expected double spin asymmetries as defined in (3.10) and 
(3.11) in the broken and unbroken (AAC) sea scenarios, as specified in fig. 1, in LO as well 
as in NLO of perturbative QCD. Here, at large scales fip ~ Mw, the perturbative stability 
with respect to sizeable variations of the factorization scale /xp is even more pronounced 
than for the Drell-Yan (7*) asymmetries in fig. 3. The expected statistical errors which 
are reduced as compared to their size in the Drell-Yan (7*) production asymmetries are 
shown in fig. 3. It is seen that in any case, LO or NLO estimates, a distinction between 
both scenarios is possible due to the sizeable reduction in the statistical errors involved 
which is mainly due to the increased luminosity at v^S = 500 GeV as compared to 
y/S = 50 GeV, relevant for Drell-Yan (7*) production as discussed in Section 3.1. We 
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recall here the envisaged integrated luminosities at RHIC, i.e. C(y/S) = (\^S/500 GeV) 
800 pb~ l and eq. (3.5). In fig. 7 we present the corresponding single asymmetries 
of eqs. (3.8) and (3.9) which turn out be less sensitive to the flavor-broken sea densities 
than the double spin asymmetries at y > in fig. 6. An interesting feature demonstrated 
here is the quality of the 'crude' approximations in eqs. (3.8) and (3.9) which are well 
satisfied at \y\ £ \ in the LO calculations where they are relevant. The ratios of singly 
and doubly polarized cross sections in eqs. (3.14) and (3.15), together with their limiting 
LO 'crude' approximations, are shown in fig. 8. Again, on account of the rather small 
expected statistical errors, it is likely that future measurements can discriminate between 
the flavor-broken (solid curves) and unbroken (dashed curves) light-sea scenarios. 

Finally we present the following ratio of combinations of polarized and unpolarized 



cross sections [26fl, which turn out to be less sensitive to absolute normalization uncer- 
tainties, 



Ao™ + + AaW- - ( Aa™ + + A^i 

pp 1 pp \ pn 1 pn I 



a J* ~ a w+ + a w- + a w+ +a w~ ' ( 3 - 17 ) 

pp PP pn 1 pn 

lo [Au{l) - Ad{l)\ [u(2) - d{2)} - [Au{l) - Ad{l)} \u{2) - d{2)} 
HI) + d(l)] [u(2) + d{2)} + [u(l) + 1(1)] K2) + d{2)} 

and its corresponding doubly polarized counterpart 
Aa^ + + Acr5T - ( AaZ + + Aaf " 



w _ P_P_ PP y pn pn / 

U pN ~ a W+ + a W~ + a W+ + a W~ y°- L °) 

pp ' pp ' pn ' " pn 

lo [Au{l) - Ad(l)\ [Au{2) - Ad{2)] - [Au{l) - Ad{\)\ [Au{2) - Ad(2)} 
[u{\) + d{l)] [u{2) + d{2)} + [u{\) + d{\)] [u{2) + d{2)} ' 

Both asymmetries have no polarized gluon contribution proportional to Ag in NLO, which 
cancels in the various differences Aa^ p — Aa^ n and Acr^- — Aa^^. Furthermore, the double 
helicity asymmetry in (3.18) vanishes for Am = Ad in LO as well as in NLO and is thus 
an interesting combination to observe the effects of a flavor-broken sea due to Au ^ Ad. 
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This is explicitly demonstrated in fig. 9 where the sizeably different expectations of the 
flavor-broken sea scenario should be easily discernible experimentally. 

Finally, for Z° production we have found that the following two double and single spin 
asymmetries are best suited for the investigation of the flavor structure of the sea: 

aZ° _ A Z° I. 

A pp = A(T PP 1 0-J 

lo a u [Au(l)Au(2) + Au(l)A«(2)] + a d \Ad(\)Ad(2) + Ad(l)Ad(2)] 
a u [u(l)u(2) + u(l)u(2)} + a d [d(l)d(2) + d(l)d(2)} 

A zo = ^pK-^4 

pp ' pn 

lo -[f3 u Au(l) - p d Ad(l)} [u(2) - d(2)} - \/3 u Au(l) - p d Ad(l)} [u(2) - d(2)} 
[a u u(l) + a d d(l)} [u{2) + d{2)] + [a u u(l) + a d d(l)} [u(2) + d{2)] 

where a q = v 2 q + a 2 q and (3 q = 2v q a q as explained and given in the Appendix. These asym- 
metries are depicted, together with their expected statistical errors, in fig. 10. It should 
be emphasized that these Z°-production asymmetries are up to an order of magnitude 
larger in the flavor-broken (Am 7^ Ad) scenario than for flavor-symmetric sea densities 
(Am = Ad) where they become almost unmeasurably small. 



(3.19) 



(3.20) 



4 Summary and Conclusions 

The possibility to determine the flavor structure of the polarized antiquark (sea) distribu- 
tions of the nucleon via vector boson (7*, W ± , Z°) production at high energy polarized 

— # 

hadron-hadron (pp, pn(d)) colliders was investigated. The perturbative stability of the 
expected asymmetries for two representative models for the flavor structure of the sea dis- 
tributions was studied and has shown that the predicted distinctive signatures for both 
flavor-symmetric and flavor-asymmetric models remain essentially unchanged in LO and 
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NLO of perturbative QCD. This demonstrates that these characteristic and distinctive 
features are genuine signatures of the models under consideration for the flavor structure 
of the polarized sea. 

In particular the polarized Drell-Yan (7*) dilepton production asymmetry AZ_ H in 
(3.2), as obtained from pp and pn collisions, or the ratio R p d of polarized pp and pn 
production cross sections in (3.3) provide us with characteristic and direct signatures for a 
flavor-broken polarized sea, Au 7^ Ad, as illustrated in fig. 3. At much larger factorization 
scales fj? F ~ M^ z , the double spin asymmetries for W ± production in (3.10) and 
(3.11), or combinations of doubly polarized pp and pn W^-production cross sections in 
(3.18) constitute similar clean and distinctive observables for studying the flavor structure 
of the polarized light sea, as shown in figs. 6 and 9. Somewhat less sensitive signatures for 
a flavor-broken sea are provided by the single spin asymmetries A^ p in (3.8) and (3.9), 
although they may give access directly to Ad(x, M^) and Au(x, M^) in specific kinematic 
regions (y <, —0.5). Such direct signatures for Ad and Au could also be obtained from 
studying ratios of W + and W~ cross sections of singly and doubly polarized pd collisions, 
cf. (3.14) and (3.16). Finally, the double spin asymmetry A— production at RHIC is an 
equally useful observable for delineating the flavor structure of the polarized sea, since it 
is expected to be about an order of magnitude larger for a flavor-broken (Au 7^ Ad) than 
for a flavor-symmetric (Au = Ad) polarized sea scenario as shown in fig. 10. 

The resolution power of the asymmetries studied depends of course on the expected 
statistical errors which were estimated for the envisaged beam polarizations and luminosi- 
ties at RHIC. They point towards the superiority of vector boson (W ± , Z°) production 
(figs. 6 - 10) over the common Drell-Yan 7* production (fig. 3) as a tool for studying the 
flavor structure of the polarized sea distributions. This derives mainly from the increased 
luminosity at the corresponding higher energies involved in vector boson production, i.e. 

~ 500 GeV, as compared to y/S ~ 50 GeV relevant for 7* (dilepton) production as 
discussed in Section 3.1. 
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Appendix 



Here we summarize all those unpolarized and polarized cross sections for (Drell-Yan) 
vector boson {i*,W ± ,Z°) production in LO and NLO(MS) needed for calculating the 
various spin-asymmetries suggested and studied in this paper. 



In terms of cross sections of definite positive and negative hadron helicities (±), an 

unpolarized cross section is generally defined by a = \{o~ ++ + o" + _ + o"_ + + a ). The 

relevant differential unpolarized Drell-Yan cross section for h\h 2 — > 7*X — > l + l~X can be 
written as 



&NP Ax F _■£->_ ' L, J„, 



q=u,d,s 



Dq^(xi, X2, Xi, x 2 ) + 2^^Q$ ^Xli X2, x±, x 2 , ^2 
x {q(xi, fip)q(x2, l4) + q(x u fj 2 F )q(x 2 , ftp)} 



+ £ jD ^ ) [ Xl > X2 > x i> x °>~2- ) 9(xi,Vf) {q(x2,^ 2 F ) + q(x 2 ,^ 2 F )} 



M 5 



„0 Jd 



M' 

2--' l!l V'' 1 "'' 2 ' 



with iV 7 * = 47ra; 2 /9S', a s = a s (/j, 2 F ), D^(xi, x 2 , X?, = — rc5)<5(rc 2 — a^O/fai + x^) 
and according to the NLO(MS) results of [28,23] 



(A.l) 



M 2 

9 1 ? \ 1' 2 ' 1' 2' 2 



T I ™0 

^ i ■ ' 2 



7T 



- 8 + 2Li 2 K) + 2U 2 {x° 2 ) 



+ ln 2 (l - x°) + ln 2 (l - x° 2 ) + 2 In In 



1 - x? 1 - x° 2 
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+ 



5(xi - x<j) 

t*0 _|_ ™0 



-02,2 
2 ~~ r~ ^ 



2 In 2 ' 2 



2^2 ^2 x\{x,2 — X 2 ) ^2 



x^ + x 2 , /ln(l -x^/x^A xS^ + x 2 , 1 
x 2 - x° J 



.X n 



+ 



x 2 , (x 2 — 



In 



(sg + sijXl-Si) 
x?(x5 + x 2 ) 



+ (!<-> 2) 



G A (xi, X 2 , Xi, X 2 ) jjAi 0\ 

+ 77 — TTT^ + H fa, X 2 , X x , X 2 ) 



[(Xt - X1)(X2 ~ X°2)} + 



+ ln 



M 2 J <5(xi - x?) <5(x 2 - s°) 



/4 



TO I ,-y»0 

1 "i - ^ 2 



3 + 2 In ^— + 2 In 



1 -xS 



XV 



.X 9 



5(xi - x?) x° 2 2 + x 2 . 1 
I -Sr — if 2 - 2 7 sr- + (l <-> 2) 

x^ + x^ x^ (x 2 -x^) 



(A.2) 



-^g^ (^1)^2)^1)^2) ,,2 



M 2 

/4 



TV 



gxg - xl) 
(x° + x^x^ 



X 



f + (xi - x?) 2 ) In 



(x° + x^)(l -x°)(xi -x°) 
x5x°(xi + x 2 ) 



+2x°(xx -x°) 



+ 



G C {xi,X 2 ,x\,X 2 

(x 2 - x%)+ 



+ H c (x 1 ,x 2 ,x { i,x {J 2 ) 



„,0\ 



+ ln^ 



M 2 j 5(x 2 - x°) h02 , , 

A*F I (^l + a '2) a 'l 



(xf + (xi - x\f 



(A.3) 



where C F = 4/3, T F = 1/2, Li 2 (x) = - Jq dthi(l - t)/t and 

(xi + x 2 )(x° 2 x° 2 + x 2 :r 2 ) 



G A (x u x 2 ,x\,x%) = 
H A (x 1 ,x 2 ,x° 1 ,x° 2 ) = 



x\x 2 2 (x\ + x 2 )(xi + x° 2 ) 
2 



XiX 2 (xi + x 2 ) 
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G c (x 1 ,x 2 ,x° 1 ,x° 2 ) 
H c (x 1 ,x 2 ,x° 1 ,x° 2 ) 



xfx 2 (a;° + x 2 ) 

Xi(x5 + x 2 )(x 2 - x°) + 2:r ( j ) :r 2 (£i + x 2 ) 
xfxl(xi + X2) 2 



(A.4) 



Furthermore 



M*\ 

X 2 , Xi, X 2 , X 1 , „ I 

/4/ 



(A.5) 



with the factorization scale /i F , usually assumed to be given by M, and x\ 2 given by (3.4), 
with the constraints r < x\ 2 < 1 and — 1 + t < < 1 — r. Alternatively one may define 
x i2 = V^ e±y with the lepton pair rapidity y in the hadron-hadron cm. system being 
constrained by In \fr < y < — In ^/t- The unpolarized cross sections for vector boson 
(1/5 = W ± , Z°) production can be written as (for simplicity we use the same symbols for 
the coefficient functions as in (A.l) although they obviously differ for different processes) 



d<l(?/,/4) = N vBy c vB [\ Xi f\ X2 
dy ^ J x o J x o 



n(°)fV „0\ , ^± / _o it-J yg 

qq v 1' 2? 2/ -^gg I J-i, x 2 , x 2 , 2 

Z7T y /i^i 



M, 2 



x {?(xi, H 2 F )q\x 2 , n 2 F ) + q(x 1 ,/i 2 F )q'(x 2 , /4)j 



+^- D £ ) ( xi,x 2 ,x?,^, ) #(xi,/4) {q'(x 2 ,H 2 F ) + q'(x 2 ,/j, F )} 



Ml 



+ 7^ D $ {^uX2,xl,x%,^^j {q(x!, fi 2 F ) + q(x u fi 2 F )} g(x 2 , fi 2 F ) 



(A.6) 



with N VB = V2nG F M^ B /3S, using M w = 80.42 GeV and M z = 91.19 GeV, and 



-qq' 



V qq '\ 2 with the relevant CKM matrix elements V U( i ~ 0.97 and V^ s ~ 0.22, and 
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c qq> = <W ( v q + a l) with v\ + S3 0.29 and ^ + w 0.37. Furthermore we have now 



Dfhxi, x 2 , x®, x°) = 5(a;i — x ( j ) )5(a;2 — ^2) an d the NLO(MS) coefficients and -Dgg j are 



(i) 



given by eqs. (15) and (16) of ref. [29], with D$ being again given by (A. 5) 



The doubly longitudinally polarized cross sections are generally defined via Ao^ ^ = 

\{o~ ++ — c + _ — o~_ + + a ). The relevant differential polarized Drell-Yan cross section for 

hih 2 — > 7*X — > l + l~X is given by (1,2 denote the arguments Xi^^f and . . . denotes the 
variables of the coefficient functions in (A.2-A.6) ) 

dAaZ\ (x F ,M 2 ,ul) ^ r 1 f 1 

dM*dx F 9 L, 7,0 



q=u,d,s 



X • 



x 



{Ag(l)Ag(2) + Ag(l)Ag(2)} 



+ ^A J D«(...)A^(l){Ag(2) + Ag(2)} 



2tt 



+ ^A J D«(...){Ag(l) + Ag(l)}A^(2) 



2tt 



where -D^' 1 ^ are as in (A.l) with given in (A. 2) and [23] 



(A.7) 



^DqJ ( Xi, X2, X 1: X 2 , 2 

A*F 



£(:r 2 - z|) 



(2x? - xi) In 



(x° + x°)(l -x°)(xi -z?) 



+2(xi - x?) 



1^2 — ^2) + 



+ ln — i / 0/ 2 ( 2x i - x 

A*F I V 3 '! X 2) X 1 



(A.8) 
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with H c '(. . . ) given in (A.4) and 

AG c (x 1 ,x 2 ,x° 1 ,x° 2 ) = 
Similarly to (A. 5) we have 



(/y»0 ,yjO \ 2 I ( ,->•» ,->•» , y*0 ,y»0 \ 2 

12/ ~~ \**- , l**-'2 12/ 

xf^Ka;] 1 + x 2 ) 



(A.9) 



AD« ^i,x 2 ,x?,^,^ = AD« ^ 2 ,x 1 ,4,x°, 



/4/ 



(A.10) 



The doubly longitudinally polarized cross sections for vector boson production are given 
by 



dy 



1,1 



X • 



)<°>( ) + 2!L D W 

'99 v- • • ) ~ 2 n 11 



x{Ag(l)Ag'(2) + Ag(l)Ag'(2)} 



+^AD^(...)Ag(l){A q '(2) + Aq'(2)} 



+ ^AD«(...){Ag(l) + Ag(l)}A^(2) 



(A.11) 



with N VB , Cq ? and -D^' 1 ^ as in (A. 6) and AD^J is now given by eq. (18) of ref. 
which relates to AD^J again via (A. 10). 



Finally, the singly longitudinally polarized cross sections for vector boson (VB = W ± , 

Z°) production are generally defined via Aa% h = \{o~ ++ + c + _ — o"_ + — a ) which are 

given by [29] 

dAaY B (y,fx 2 F ) ^ ri r 1 



dy 



q,q' 
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X • 



D 



(0), 
qq \ 



2tt 



qq 



x{-Ag(l)g'(2) + Ag(l)g'(2)} 



+ S^ ) (---){-Ag(l)+Ag(l)}^(2) 



2tt 



(A.12) 



where all normalizations, couplings and coefficient functions are as in (A. 11) except for 
Cg°, which is now given by c^°, = 5 qq > 2v q a q , i.e. cf ° = | — | sin 2 6^ 0.19 and cf d ° = 
cf = I_| sin 2 6^ « 0.35. 
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Figure Captions 



Fig. 1. Polarized (anti)quark distributions in LO and NLO according to eqs. (2.2) - (2.6) in 
the flavor-broken sea scenario (full curves) at various fixed factorization scales up = 
Q 2 . For comparison the corresponding 'standard' flavor-symmetric sea distribution 
Aq according to AAC [ITS] is shown by the dashed curves. The remaining Am and 
Ad densities are similar in both scenarios (cf. fig. 5). 

Fig. 2. The unpolarized Drell-Yan dilepton production cross section d 2 o1 v ^ M ^ j dM dxp for 
= 50 and 100 GeV in NLO based on the GRV98 El NLO parton distributions. 



Fig. 3. The ratio R p + n of polarized pp and pn Drell-Yan cross sections in eq. (3.3) and the 
asymmetry in eq. (3.2) in LO and NLO based on the distributions of the 

broken and unbroken (AAC) sea scenarios as specified in fig. 1. The perturbative 
stability of the predicted asymmetries and their sensitivity to the choice of the 
factorization scale fip is examined in the two lower figures. Statistical errors are 
calculated according to eqs. (3.5) and (3.6). The 'crude approx.' curve refers to 
[1 + Ad(2)/A«(2)]/2 in eq. (3.3). 

Fig. 4. The unpolarized differential cross sections da™ > z °/dy for VS = 500 GeV (RHIC) 



in NLO based on the GRV98 [21] NLO parton distributions 



Fig. 5. The LO (anti)quark asymmetries (Aq/q) Aq/q at the relevant factorization scale 
H F = Q 2 = M$y, dominating the single spin asymmetries for AW p at \y\ > | 
according to eqs. (3.8) and (3.9). 

Fig. 6. The double spin asymmetries 

[eq. (3.10)] and Aj ff [eq. (3.11)] at ^/S = 500 
GeV in the broken and unbroken (AAC) sea scenarios. The statistical errors are 
evaluated according to eq. (3.5). The perturbative stability of the predicted asym- 
metries and their sensitivity to the choice of the factorization scale fip is shown as 
well. 
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Fig. 7. The single spin asymmetries Aj p + [eq. (3.8)] and A\Y p [eq. (3.9)] at V~S = 500 GeV 
in the broken and unbroken (AAC) sea scenarios, as specified in fig. 1. The sta- 
tistical errors are evaluated according to eq. (3.7). The errors on the AAC curves 
are similarly small as in the broken scenario. The quality of the 'crude' approxima- 
tions in eq. (3.8), AJ(1)/J(1), and eq. (3.9), A-u(l)/w(l), is examined as well. The 
perturbative stability of the predicted asymmetries is similar to the one shown in 
fig. 6. 

Fig. 8. The ratios Aa^ + /Acx~ and Acr^ + / Ao^ d for doubly and singly polarized cross 
sections, respectively, at = 500 GeV obtained via eqs. (3.15) and (3.14). The 
meaning of the curves corresponds to that in fig. 7. 

Fig. 9. The asymmetries and o^ N at = 500 GeV in eqs. (3.18) and (3.17), respec- 
tively, for the broken and unbroken (AAC) scenario. 

Fig. 10. The asymmetries Af- and Af p _ n in LO and NLO at = 500 GeV given in eqs. 

(3.19) and (3.20), respectively. The meaning of the curves corresponds to that in 
fig. 9. 
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